The exon junction complex is required for DMD gene splicing fidelity and myogenic differentiation.

Deposition of the exon junction complex (EJC) upstream of exon-exon junctions helps maintain transcriptome integrity by preventing spurious re-splicing events in already spliced mRNAs. Here we investigate the importance of EJC for the correct splicing of the 2.2-megabase-long human DMD pre-mRNA, which encodes dystrophin, an essential protein involved in cytoskeletal organization and cell signaling. Using targeted RNA-seq, we show that knock-down of the eIF4A3 and Y14 core components of EJC in a human muscle cell line causes an accumulation of mis-splicing events clustered towards the 3' end of the DMD transcript (Dp427m). This deregulation is conserved in the short Dp71 isoform expressed ubiquitously except in adult skeletal muscle and is rescued with wild-type eIF4A3 and Y14 proteins but not with an EJC assembly-defective mutant eIF4A3. MLN51 protein and EJC-associated ASAP/PSAP complexes independently modulate the inclusion of the regulated exons 71 and 78. Our data confirm the protective role of EJC in maintaining splicing fidelity, which in the DMD gene is necessary to preserve the function of the critical C-terminal protein-protein interaction domain of dystrophin present in all tissue-specific isoforms. Given the role of the EJC in maintaining the integrity of dystrophin, we asked whether the EJC could also be involved in the regulation of a mechanism as complex as skeletal muscle differentiation. We found that eIF4A3 knockdown impairs myogenic differentiation by blocking myotube formation. Collectively, our data provide new insights into the functional roles of EJC in human skeletal muscle.

RFC1 nonsense and frameshift variants cause CANVAS: clues for an unsolved pathophysiology.

Cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) is an inherited late-onset neurological disease caused by bi-allelic AAGGG pentanucleotide expansions within intron 2 of RFC1. Despite extensive studies, the pathophysiological mechanism of these intronic expansions remains elusive. We screened by clinical exome sequencing two unrelated patients presenting with late-onset ataxia. A repeat-primer polymerase chain reaction was used for RFC1 AAGGG intronic expansion identification. RFC1 mRNA expression was assessed by quantitative reverse transcription-polymerase chain reaction. We identified the first two CANVAS affected patients who are compound heterozygous for RFC1 truncating variants (p.Arg388* and c.575delA, respectively) and a pathological AAGGG expansion. RFC1 expression studies in whole blood showed a significant reduction of RFC1 mRNA for both patients compared to three patients with bi-allelic RFC1 expansions. In conclusion, this observation provides clues that suggest bi-allelic RFC1 conditional loss-of-function as the cause of the disease.

Association between prophylactic angiotensin-converting enzyme inhibitors and overall survival in Duchenne muscular dystrophy-analysis of registry data.

AIMS: To estimate the effect of prophylactic angiotensin-converting enzyme inhibitors (ACEi) on survival in Duchenne muscular dystrophy (DMD). METHODS AND RESULTS: We analysed the data from the French multicentre DMD Heart Registry (ClinicalTrials.gov: NCT03443115). We estimated the association between the prophylactic prescription of ACEi and event-free survival in 668 patients aged 8 to 13 years, with normal left ventricular function, using (i) a Cox model with intervention as a time-dependent covariate, (ii) a propensity-based analysis comparing ACEi treatment vs. no treatment, and (iii) a set of sensitivity analyses. The study outcomes were overall survival and hospitalizations for heart failure (HF) or acute respiratory failure. Among the 668 patients included in the DMD Heart Registry, 576 (mean age 6.1 ± 2.8 years) were eligible for this study, of whom 390 were treated with ACEi prophylactically. Death occurred in 53 patients (13.5%) who were and 60 patients (32.3%) who were not treated prophylactically with ACEi, respectively. In a Cox model with intervention as a time-dependent variable, the hazard ratio (HR) associated with ACEi treatment was 0.49 [95% confidence interval (CI) 0.34-0.72] and 0.47 (95% CI 0.31-0.17) for overall mortality after adjustment for baseline variables. In the propensity-based analysis, 278 patients were included in the treatment group and 834 in the control group, with 18.5% and 30.4% 12-year estimated probability of death, respectively. ACEi were associated with a lower risk of death (HR 0.39; 95% CI 0.17-0.92) and hospitalization for HF (HR 0.16; 95% CI 0.04-0.62). All other sensitivity analyses yielded similar results. CONCLUSION: Prophylactic ACEi treatment in DMD was associated with a significantly higher overall survival and lower rates of hospitalization for HF.

First Identification of RNA-Binding Proteins That Regulate Alternative Exons in the Dystrophin Gene.

The Duchenne muscular dystrophy (DMD) gene has a complex expression pattern regulated by multiple tissue-specific promoters and by alternative splicing (AS) of the resulting transcripts. Here, we used an RNAi-based approach coupled with DMD-targeted RNA-seq to identify RNA-binding proteins (RBPs) that regulate splicing of its skeletal muscle isoform (Dp427m) in a human muscular cell line. A total of 16 RBPs comprising the major regulators of muscle-specific splicing events were tested. We show that distinct combinations of RBPs maintain the correct inclusion in the Dp427m of exons that undergo spatio-temporal AS in other dystrophin isoforms. In particular, our findings revealed the complex networks of RBPs contributing to the splicing of the two short DMD exons 71 and 78, the inclusion of exon 78 in the adult Dp427m isoform being crucial for muscle function. Among the RBPs tested, QKI and DDX5/DDX17 proteins are important determinants of DMD exon inclusion. This is the first large-scale study to determine which RBP proteins act on the physiological splicing of the DMD gene. Our data shed light on molecular mechanisms contributing to the expression of the different dystrophin isoforms, which could be influenced by a change in the function or expression level of the identified RBPs.

Tumor Necrosis Factor Receptor SF10A (TNFRSF10A) SNPs Correlate With Corticosteroid Response in Duchenne Muscular Dystrophy.

BACKGROUND: Duchenne muscular dystrophy (DMD) is a rare and severe X-linked muscular dystrophy in which the standard of care with variable outcome, also due to different drug response, is chronic off-label treatment with corticosteroids (CS). In order to search for SNP biomarkers for corticosteroid responsiveness, we genotyped variants across 205 DMD-related genes in patients with differential response to steroid treatment. METHODS AND FINDINGS: We enrolled a total of 228 DMD patients with identified dystrophin mutations, 78 of these patients have been under corticosteroid treatment for at least 5 years. DMD patients were defined as high responders (HR) if they had maintained the ability to walk after 15 years of age and low responders (LR) for those who had lost ambulation before the age of 10 despite corticosteroid therapy. Based on interactome mapping, we prioritized 205 genes and sequenced them in 21 DMD patients (discovery cohort or DiC = 21). We identified 43 SNPs that discriminate between HR and LR. Discriminant Analysis of Principal Components (DAPC) prioritized 2 response-associated SNPs in the TNFRSF10A gene. Validation of this genotype was done in two additional larger cohorts composed of 46 DMD patients on corticosteroid therapy (validation cohorts or VaC1), and 150 non ambulant DMD patients and never treated with corticosteroids (VaC2). SNP analysis in all validation cohorts (N = 207) showed that the CT haplotype is significantly associated with HR DMDs confirming the discovery results. CONCLUSION: We have shown that TNFRSF10A CT haplotype correlates with corticosteroid response in DMD patients and propose it as an exploratory CS response biomarker.

Clinical Phenotypes of DMD Exon 51 Skip Equivalent Deletions: A Systematic Review.

BACKGROUND: Eteplirsen, the first FDA-approved RNA-modifying therapy for DMD, is applicable to ∼13% of patients with DMD. Because multiple exonic deletions are amenable to exon 51 skipping, the isoforms resulting from the various exon 51-skipped transcripts may vary in stability, function, and phenotype. OBJECTIVE/METHODS: We conducted a detailed review of dystrophinopathy published literature and unpublished databases to compile phenotypic features of patients with exon 51 "skip-equivalent" deletions. RESULTS: Theoretically, 48 different in-frame transcripts may result from exon 51 skipping. We found sufficient clinical information on 135 patients carrying mutations that would result in production of 11 (23%) of these transcripts, suggesting the remainder have not been identified in vivo. The majority had mild phenotypes: BMD (n = 81) or isolated dilated cardiomyopathy (n = 3). Particularly interesting are the asymptomatic (n = 10) or isolated hyperCKemia (n = 20) patients with deletions of exons 45- 51, 48- 51, 49- 51 and 50- 51. Finally, 16 (12%) had more severe phenotypes described as intermediate (n = 2) or DMD (n = 14), and 6 reports had no definitive phenotype. CONCLUSIONS: This review shows that the majority of exon 51 "skip-equivalent" deletions result in milder (BMD) phenotypes and supports that exon 51 skipping therapy could provide clinical benefit, although we acknowledge that other factors, such as age at treatment initiation or ongoing standard of care, may influence the degree of benefit.

EMQN best practice guidelines for genetic testing in dystrophinopathies.

Dystrophinopathies are X-linked diseases, including Duchenne muscular dystrophy and Becker muscular dystrophy, due to DMD gene variants. In recent years, the application of new genetic technologies and the availability of new personalised drugs have influenced diagnostic genetic testing for dystrophinopathies. Therefore, these European best practice guidelines for genetic testing in dystrophinopathies have been produced to update previous guidelines published in 2010.These guidelines summarise current recommended technologies and methodologies for analysis of the DMD gene, including testing for deletions and duplications of one or more exons, small variant detection and RNA analysis. Genetic testing strategies for diagnosis, carrier testing and prenatal diagnosis (including non-invasive prenatal diagnosis) are then outlined. Guidelines for sequence variant annotation and interpretation are provided, followed by recommendations for reporting results of all categories of testing. Finally, atypical findings (such as non-contiguous deletions and dual DMD variants), implications for personalised medicine and clinical trials and incidental findings (identification of DMD gene variants in patients where a clinical diagnosis of dystrophinopathy has not been considered or suspected) are discussed.

TCTEX1D1 is a genetic modifier of disease progression in Duchenne muscular dystrophy.

Duchenne muscular dystrophy (DMD) is caused by pathogenic variants in the DMD gene leading to the lack of dystrophin. Variability in the disease course suggests that other factors influence disease progression. With this study we aimed to identify genetic factors that may account for some of the variability in the clinical presentation. We compared whole-exome sequencing (WES) data in 27 DMD patients with extreme phenotypes to identify candidate variants that could affect disease progression. Validation of the candidate SNPs was performed in two independent cohorts including 301 (BIO-NMD cohort) and 109 (CINRG cohort of European ancestry) DMD patients, respectively. Variants in the Tctex1 domain containing 1 (TCTEX1D1) gene on chromosome 1 were associated with age of ambulation loss. The minor alleles of two independent variants, known to affect TCTEX1D1 coding sequence and induce skipping of its exon 4, were associated with earlier loss of ambulation. Our data show that disease progression of DMD is affected by a new locus on chromosome 1 and demonstrate the possibility to identify genetic modifiers in rare diseases by studying WES data in patients with extreme phenotypes followed by multiple layers of validation.

Current management of Duchenne muscular dystrophy in the Middle East: expert report.

Aim: Duchenne muscular dystrophy (DMD) is a severe and rare X-linked neuromuscular childhood disorder that results in functional decline, loss of ambulation and early death due to cardiac or respiratory failure. The objective of this paper is to address different aspects of the current management of DMD in the Middle East, north Africa (MENA) region, and to gather experts' recommendations on how to optimally diagnose and treat patients suffering from this disease. Methods: A group of experts (neuromuscular medicine, neuropediatricians and geneticists) convened to discuss the diagnosis and management of DMD in the MENA region. A list of practical statements was prepared by the chair of the meeting to guide the discussions around critical aspects relating to the current and future management of DMD. Results & conclusion: Ideally, DMD management should be a multidisciplinary approach. Nevertheless, few tertiary care hospitals in the region are currently able to provide the full spectrum of medical expertise and services needed by DMD patients. Clinical practice in the region remains heterogeneous. Specific guidelines for diagnosis and treatment are needed in the MENA region to improve outcomes. Disease awareness among the general public and the medical community is lacking. Now that mutation-specific therapies are being developed and more widely studied, general education programs regarding early signs and symptoms, a standardized referral and diagnosis pathway, patient registries and support groups will significantly improve the management of the disease.

MoBiDiC Prioritization Algorithm, a Free, Accessible, and Efficient Pipeline for Single-Nucleotide Variant Annotation and Prioritization for Next-Generation Sequencing Routine Molecular Diagnosis.

Interpretation of next-generation sequencing constitutes the main limitation of molecular diagnostics. In diagnosing myopathies and muscular dystrophies, another issue is efficiency in predicting the pathogenicity of variants identified in large genes, especially TTN; current in silico prediction tools show limitations in predicting and ranking the numerous variants of such genes. We propose a variant-prioritization tool, the MoBiDiCprioritization algorithm (MPA). MPA is based on curated interpretation of data on previously reported variants, biological assumptions, and splice and missense predictors, and is used to prioritize all types of single-nucleotide variants. MPA was validated by comparing its sensitivity and specificity to those of dbNSFP database prediction tools, using a data set composed of DYSF, DMD, LMNA, NEB, and TTN variants extracted from expert-reviewed and ExAC databases. MPA obtained the best annotation rates for missense and splice variants. As MPA aggregates the results from several predictors, individual predictor errors are counterweighted, improving the sensitivity and specificity of missense and splice variant predictions. We propose a sequential use of MPA, beginning with the selection of variants with higher scores and followed by, in the absence of candidate pathologic variants, consideration of variants with lower scores. We provide scripts and documentation for free academic use and a validated annotation pipeline scaled for panel and exome sequencing to prioritize single-nucleotide variants from a VCF file.

Identification of Splicing Factors Involved in DMD Exon Skipping Events Using an In Vitro RNA Binding Assay.

Mutation-induced exon skipping in the DMD gene can modulate the severity of the phenotype in patients with Duchenne or Becker Muscular Dystrophy. These alternative splicing events are most likely the result of changes in recruitment of splicing factors at cis-acting elements in the mutated DMD pre-mRNA. The identification of proteins involved can be achieved by an affinity purification procedure. Here, we provide a detailed protocol for the in vitro RNA binding assay that we routinely apply to explore molecular mechanisms underlying splicing defects in the DMD gene. In vitro transcribed RNA probes containing either the wild type or mutated sequence are oxidized and bound to adipic acid dihydrazide-agarose beads. Incubation with a nuclear extract allows the binding of nuclear proteins to the RNA probes. The unbound proteins are washed off and then the specifically RNA-bound proteins are released from the beads by an RNase treatment. After separation by SDS-PAGE, proteins that display differential binding affinities for the wild type and mutant RNA probes are identified by mass spectrometry.

Normal and altered pre-mRNA processing in the DMD gene.

Splicing of pre-mRNA is a crucial regulatory stage in the pathway of gene expression controlled by multiple post- and co-transcriptional mechanisms. The large Duchenne muscular dystrophy gene encoding the protein dystrophin provides a striking example of the complexity of human pre-mRNAs. In this review, we summarize the current state of knowledge about canonical and non-canonical splicing in the DMD pre-mRNA, with a focus on mechanisms that take place in the full-length transcript isoform expressed in human skeletal muscle. In particular, we highlight recent work demonstrating that multi-step events are required for long DMD intron removal. The role of temporary intron retention in the occurrence of alternative splicing events is also discussed. Even though the proportion of splicing mutations is lower than reported in other genes, a great diversity of splicing defects linked to point mutations, but also large genomic rearrangements are observed in the DMD gene. We provide an overview of the molecular mechanisms underlying aberrant splicing in patients with Duchenne or Becker muscular dystrophy, and we also detail how alternative splicing can serve as a disease modifier in patients by changing the outcome of the primary defect.

Targeted RNA-Seq profiling of splicing pattern in the DMD gene: exons are mostly constitutively spliced in human skeletal muscle.

We have analysed the splicing pattern of the human Duchenne Muscular Dystrophy (DMD) transcript in normal skeletal muscle. To achieve depth of coverage required for the analysis of this lowly expressed gene in muscle, we designed a targeted RNA-Seq procedure that combines amplification of the full-length 11.3 kb DMD cDNA sequence and 454 sequencing technology. A high and uniform coverage of the cDNA sequence was obtained that allowed to draw up a reliable inventory of the physiological alternative splicing events in the muscular DMD transcript. In contrast to previous assumptions, we evidenced that most of the 79 DMD exons are constitutively spliced in skeletal muscle. Only a limited number of 12 alternative splicing events were identified, all present at a very low level. These include previously known exon skipping events but also newly described pseudoexon inclusions and alternative 3' splice sites, of which one is the first functional NAGNAG splice site reported in the DMD gene. This study provides the first RNA-Seq-based reference of DMD splicing pattern in skeletal muscle and reports on an experimental procedure well suited to detect condition-specific differences in this low abundance transcript that may prove useful for diagnostic, research or RNA-based therapeutic applications.

Association Study of Exon Variants in the NF-κB and TGFß Pathways Identifies CD40 as a Modifier of Duchenne Muscular Dystrophy.

The expressivity of Mendelian diseases can be influenced by factors independent from the pathogenic mutation: in Duchenne muscular dystrophy (DMD), for instance, age at loss of ambulation (LoA) varies between individuals whose DMD mutations all abolish dystrophin expression. This suggests the existence of trans-acting variants in modifier genes. Common single nucleotide polymorphisms (SNPs) in candidate genes (SPP1, encoding osteopontin, and LTBP4, encoding latent transforming growth factor ß [TGFß]-binding protein 4) have been established as DMD modifiers. We performed a genome-wide association study of age at LoA in a sub-cohort of European or European American ancestry (n = 109) from the Cooperative International Research Group Duchenne Natural History Study (CINRG-DNHS). We focused on protein-altering variants (Exome Chip) and included glucocorticoid treatment as a covariate. As expected, due to the small population size, no SNPs displayed an exome-wide significant p value (< 1.8 × 10-6). Subsequently, we prioritized 438 SNPs in the vicinities of 384 genes implicated in DMD-related pathways, i.e., the nuclear-factor-κB and TGFß pathways. The minor allele at rs1883832, in the 5'-untranslated region of CD40, was associated with earlier LoA (p = 3.5 × 10-5). This allele diminishes the expression of CD40, a co-stimulatory molecule for T cell polarization. We validated this association in multiple independent DMD cohorts (United Dystrophinopathy Project, Bio-NMD, and Padova, total n = 660), establishing this locus as a DMD modifier. This finding points to cell-mediated immunity as a relevant pathogenetic mechanism and potential therapeutic target in DMD.

A recurrent deletion in the SLC5A2 gene including the intron 7 branch site responsible for familial renal glucosuria.

Familial renal glycosuria (FRG) is caused by mutations in the SLC5A2 gene, which codes for Na+-glucose co-transporters 2 (SGLT2). The aim of this study was to analyze and identify the mutations in 16 patients from 8 families with FRG. All coding regions, including intron-exon boundaries, were analyzed using PCR followed by direct sequence analysis. Six mutations in SLC5A2 gene were identified, including five missense mutations (c.393G > C, p.K131N; c.1003A > G, p.S335G; c.1343A > G, p.Q448R; c.1420G > C, p.A474P; c.1739G > A, p.G580D) and a 22-bp deletion in intron 7 (c.886(-10_-31)del) removing the putative branch point sequence. By the minigene studies using the pSPL3 plasmids, we confirmed that the deletion c.886(-10_-31)del acts as a splicing mutation. Furthermore, we found that this deletion causes exclusion of exon 8 in the SCL5A2 transcript in patients. The mutation c.886(-10_-31)del was present in 5 (62.5%) of 8 families, and accounts for about 37.5% of the total alleles (6/16). In conclusion, six mutations resulting in FRG were found, and the c.886(-10_-31)del may be the high frequency mutation that can be screened in FRG patients with uniallelic or negative SLC5A2 mutations.

Two Novel HOGA1 Splicing Mutations Identified in a Chinese Patient with Primary Hyperoxaluria Type 3.

BACKGROUND: Twenty-six HOGA1 mutations have been reported in primary hyperoxaluria (PH) type 3 (PH3) patients with c.700 + 5G>T accounting for about 50% of the total alleles. However, PH3 has never been described in Asians. METHODS: A Chinese child with early-onset nephrolithiasis was suspected of having PH. We searched for AGXT, GRHPR and HOGA1 gene mutations in this patient and his parents. All coding regions, including intron-exon boundaries, were analyzed using PCR followed by direct sequence analysis. RESULTS: Two heterozygous mutations not previously described in the literature about HOGA1 were identified (compound heterozygous). One mutation was a successive 2 bp substitution at the last nucleotide of exon 6 and at the first nucleotide of intron 6, respectively (c.834_834 + 1GG>TT), while the other one was a guanine to adenine substitution of the last nucleotide of exon 6 (c.834G>A). Direct sequencing analysis failed to find these mutations in 100 unrelated healthy subjects and the functional role on splicing of both variants found in this study was confirmed by a minigene assay based on the pSPL3 exon trapping vector. In addition, we found a SNP in this family (c.715G>A, p.V239I). There were no mutations detected in AGXT and GRHPR. CONCLUSION: Two novel HOGA1 mutations were identified in association with PH3. This is the first description and investigation on mutant gene analysis of PH3 in an Asian.

FUBP1: a new protagonist in splicing regulation of the DMD gene.

We investigated the molecular mechanisms for in-frame skipping of DMD exon 39 caused by the nonsense c.5480T>A mutation in a patient with Becker muscular dystrophy. RNase-assisted pull down assay coupled with mass spectrometry revealed that the mutant RNA probe specifically recruits hnRNPA1, hnRNPA2/B1 and DAZAP1. Functional studies in a human myoblast cell line transfected with DMD minigenes confirmed the splicing inhibitory activity of hnRNPA1 and hnRNPA2/B1, and showed that DAZAP1, also known to activate splicing, acts negatively in the context of the mutated exon 39. Furthermore, we uncovered that recognition of endogenous DMD exon 39 in muscle cells is promoted by FUSE binding protein 1 (FUBP1), a multifunctional DNA- and RNA-binding protein whose role in splicing is largely unknown. By serial deletion and mutagenesis studies in minigenes, we delineated a functional intronic splicing enhancer (ISE) in intron 38. FUBP1 recruitment to the RNA sequence containing the ISE was established by RNA pull down and RNA EMSA, and further confirmed by RNA-ChIP on endogenous DMD pre-mRNA. This study provides new insights about the splicing regulation of DMD exon 39, highlighting the emerging role of FUBP1 in splicing and describing the first ISE for constitutive exon inclusion in the mature DMD transcript.

Small-scale high-throughput sequencing-based identification of new therapeutic tools in cystic fibrosis.

PURPOSE: Although 97-99% of CFTR mutations have been identified, great efforts must be made to detect yet-unidentified mutations. METHODS: We developed a small-scale next-generation sequencing approach for reliably and quickly scanning the entire gene, including noncoding regions, to identify new mutations. We applied this approach to 18 samples from patients suffering from cystic fibrosis (CF) in whom only one mutation had hitherto been identified. RESULTS: Using an in-house bioinformatics pipeline, we could rapidly identify a second disease-causing CFTR mutation for 16 of 18 samples. Of them, c.1680-883A>G was found in three unrelated CF patients. Analysis of minigenes and patients' transcripts showed that this mutation results in aberrantly spliced transcripts because of the inclusion of a pseudoexon. It is located only three base pairs from the c.1680-886A>G mutation (1811+1.6kbA>G), the fourth most frequent mutation in southwestern Europe. We next tested the effect of antisense oligonucleotides targeting splice sites on these two mutations on pseudoexon skipping. Oligonucleotide transfection resulted in the restoration of the full-length, in-frame CFTR transcript, demonstrating the effect of antisense oligonucleotide-induced pseudoexon skipping in CF. CONCLUSION: Our data confirm the importance of analyzing noncoding regions to find unidentified mutations, which is essential to designing targeted therapeutic approaches.